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I. INTRODUCTION 



The dynamics of spin-glasses and other disordered systems exhibits a very much studied phenomenon known as 
"aging": the behaviour of the system depends on its history, and experiments show a typical out of equilibrium 
regime on all (accessible) time scales jjj. In the simplest case one quenches the system into its low temperature 
phase at time t = 0, and the dynamics of the system depends on its age, i.e. the time elapsed since the quench. 
This type of behavior can be studied for example by looking at the correlation function of some local observable 
0(t), C(t,t') —< 0(t)0(t') >, or at the response of such an observable to a change in a conjugated external field 
h(t'): r(t,t') =< >• While in the usual equilibrium behaviour these two-times quantities obey time translational 

invariance (TTI) (C(t,t ! ) = C(t—t'), r(t,t') — r(t—t')) and fluctuation-dissipation theorem (FDT) relating correlation 
and response, one frequently observes in off equilibrium dynamics a dependence in C(t,t') ~ t~ a C(t' /t), which is 
referred to as an aging behaviour, and violation of FDT. 

This kind of aging behaviour is not restricted to spin glasses: the persistence of out of equilibrium effects even after 
very long times has been observed in many other systems, either experimental systems [Q, or in computer simulations 
H . In some cases aging could be studied analytically || |J . 

The kind of loose definition of aging that we have used so far seems to be ubiquitous and to hide a variety of very 
distinct physical situations. While the mean field spin glass is known to be a complicated system with a rough free 
energy landscape with many metastable states, aging also occurs in much simpler problems like the random walk [Q, 
the coarsening of domain walls in a ferromagnet quenched below its critical temperature JTo) ] , or some problems with 
purely entropic barriers all problems in which the free energy landscape seems to be very simple. 

It is interesting to find a way to distinguish between these different types of aging, and this paper wants to make 
some step towards such a classification. A first classification of aging has already been proposed in the litterature || . 
In the mean field dynamics of spin glasses, it has been shown that the response function exhibits an anomaly in the 
low temperature phase [l^Ji^l . While it looks mysterious in the framework of equilibrium dynamics, this anomaly is 
well understood if one studies off equilibrium dynamics Q . The anomaly is there defined as 

r-t />00 

\ = lim / dt'r(t,t')- / lim r(t w +T,t w )dr (1) 

It measures the difference between the susceptibility of the system at large times and the susceptibility of an hypo- 
thetical system which would be at equilibrium. A non-zero anomaly shows the existence of a long term memory of 
the system to some perturbations occuring at any time. Systems with such an anomaly certainly exhibit strong aging 
effects. 

In spite of its nice mathematical structure, the anomaly is in general not easy to control and compute (analytically 
or numerically) . In this paper we want to propose another tool for the classification of aging. We shall use an overlap 
Qt m (t w + i> t w + t) between two identical copies of the system, which are constrained to evolve from the same initial 
configuration and with the same thermal noise between the initial quench and a time t w , and then evolve with different 
realizations of the thermal noise between t w and t w + 1. This quantity was introduced in |T^ ], and in a study of the 
aging dynamics of the spherical spin glass by Cugliandolo and Dean |15| (slightly different objects involving two copies 
of the system evolving with the same noise had also been studied before pH ). We argue that the asymptotic value of 
this overlap in the double limit lim tw ^oali'mt^ooQt w (tw +t,t w +t) distinguishes between different types of aging. In 
a first class of systems the limit is finite (equal to the Edwards Anderson parameter qEA in the cases we have studied 
so far). This class which we call of type I includes the models with coarsening of domain walls: we show it explicitly 
hereafter in the case of the 0(n) model with n — > oo, and within some widely used assumptions for the domain 
growth in the non conserved scalar order parameter case. The second class, aging systems of type II, contains the 
spin glass like problems with complicated free energy landscapes, and we study explicitely the p spin spherical models 
or the zero dimensional version of the manifolds in random potential. For this class, the limit limt w —>oalimt^>oo of 
the overlap is equal to the minimum possible overlap (i.e. zero for the p spin spherical model -with p > 2- and go f° r 
the zero dimensional manifolds). 

Besides suggesting a first (rough) classification, this overlap function may turn to give some intuitive ideas about 
the energy landscape in which the system evolves, and its complexity. For example, if we think of a system falling 
down a "gutter" , it is clear that it will continuously go away from its position at t w (the correlation function decreases 
to zero), but two copies separated at t w will not be able to separate indefinitely, and the overlap will have a limit at 
long times which can depend on t w : as t w grows, the system is going closer and closer to equilibrium. Type I systems 
seem to have such a behaviour. 

On the contrary, a rugged landscape with many bifurcations, and many different paths, will allow two copies to 
really move away from one another, so the overlap will decay to its minimum possible value, for any finite t w : the 
distance between the two replicas becomes the largest possible one. 
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The classification induced by the asymptotic value of the overlap function might coincide with the one quoted above, 
using the anomaly of the response. Indeed all type II systems that we study are known to possess a non vanishing 
anomaly. In type I systems, the anomaly has been computed so far only in the 0(n) model with n — > oo, where it 
does vanish (or equivalently for the p = 2 spherical model p5|| ). The existence of a general relation between these 
two criterions remains to be studied. At an intuitive level it may look plausible: if one thinks of a type I system 
as evolving in a phase space gutter, it should not have a long term memory of a perturbation. On the other hand, 
a type II system evolving in a rugged landscape will be continuously bifurcating and a change of direction will be 
remembered at long times. As we are aware that such intuitive arguments can be very misleading, we just mention 
them here as a motivation to further studies of the response anomaly in various aging systems. 

The paper is organized as follows: in section II, we define the dynamics and various quantities we study, and present 
the general features of equilibrium dynamics. Section III is devoted to the study of various problems of domain-growth, 
with analytical and numerical results. Type II systems are studied in section IV, where we analyse in particular the 
behaviour of the zero dimensional version of the random manifold problem and of Bouchaud's model of phase space 
traps The last section contains our conclusions. 



II. DEFINITIONS, EQUILIBRIUM DYNAMICS 

We consider systems described by a field 4>{x) in a d dimensional space (we shall also consider spin systems, 
with obvious generalizations of the definitions). Given a Hamiltonian H[<f>] — f d x?i((f)(x)), we assume a Langevin 
dynamics at temperature T: 

dd>(x,t) dH 

where 77 is a white noise, with < rj(x,t)r](x' ,t') >= 2T5 d (x — x')8(t — t') (<> means an average over this thermal 
noise) . 

The quantities we are mostly interested in arc the following: 

• the two-time autocorrelation function C(t,t'): it is the mean overlap between the configurations of the field at 
times t and t', 

C(t,t') = y I d d x< <f>(x,t)<f>(x,l!) >; (3) 



the response function 



• the overlap function Q tvo (t,t'): the system evolves during a certain time t w ; at t w a copy is made, and the 
two systems obtained, labeled by W and evolve independently; Qt w (t,t') is then the overlap between the 
configuration of one copy at time t and the other at time t'\ 

Q tw {tJ) = ±Jd d x< 4>^(x,t)cf>^(x,t') > . (5) 

Of course, for t < t w or t' < t w , Q tw (t, t') = C(t, t'). 

Before turning to out of equilibrium dynamics, let us first show that the overlap Qt^{t,t') is simply related to the 
correlation in the case of equilibrium dynamics. 

If a system is evolving among a set of states, according to a master equation, with transitions rates obeying detailed 
balance, i.e.: 



dt 



T ijP ?=T jip e« (6) 
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where Pi(t) is the probability of being in state i at time t, the formal solution is 

f^)=]T<z|e rt |j> Pj (0), (7) 

j 

where < i\e Tt \j > are the matrix elements of the evolution operator e Tt (< i\T\j >= Ty). The detailed balance 
implies that: 

<j\e Tt \i>pf=<i\e Tt \j>pf (8) 

If we express this property in terms of the overlap between two replicas evolving in equilibrium dynamics, we obtain 
(see figure (1)): 

Qas{s,t) =Cas(s + t) (9) 

where we have defined Q as {s, t) — lim^^oo Q tw (t w + s,t w +t) and C as (t) = lim^^oo C(t w +t,t w ). 



(1) =(2) 




Q(t,s) 



FIG. 1. the equilibrium value of the overlap, when t w — > oo between the first replica at time t w + s (point B) and the second 
at time t w +t (point C) is the same as the overlap between D (time t w — s, before the separation) and C), or between E (time 
t w — t) and B 



There exist other interesting large time limits in the problem which exhibit aging. In particular, the interesting 
property of weak-ergodicity breaking prj| , defined by lim^oo C(t w +t,t w ) — 0, expresses the fact that such system 
never reach equilibrium. In the following we will therefore be interested in the function 

S(t w ) ee lim Q tw {t w + t, t w + t) (10) 

t — >oo 

and in particular in its large t w limit Soo — limt^^oo limt^oo Q tw (t w +t, t w +t). We shall show that this limit depends 
on the type of system one considers and allows for a distinction of various aging types. 



III. DOMAIN-GROWTH PROCESSES 



A phenomenon which is often considered as the typical example of an out of equilibrium dynamical evolution 
is the phase ordering kinetics (l^Jl^]. It is the domain growth process for an infinite system with different low 
temperature ordered phases, suddenly quenched from a disordered high temperature region into an unstable state at 
low temperature. Here we shall keep to the dynamical evolution of systems with a non conserved order parameter 

EE- 

The system we study is described by a n component vector field <f>(x,t) representing the density of magnetization 
at the point x of a d dimensional space, as a function of time. The system is prepared at high temperature, where 
< <j) > =0 and then rapidly quenched a t = in a low temperature region, where there are more than one 
energetically favourable states with < <f> >=/= 0. This situation is well described by a typical coarse grained free energy: 

F = Jd d x [i(V^) 2 +V(0)], (11) 

where the first term represents the energy cost of an interface between two different phases and V[<f>] is a potential 
with minima at different values of (j). The state with < <j) >= is unstable at low temperature, so the system evolves 
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by forming larger and larger domains of a single phase; at late stage of growth, the typical pattern of domains is self 
similar and the characteristic size of a domain is L(t). This evolution can be studied for instance through a Langevin 
dynamics with thermal noise r\: 

^ = V 2 0-V'(0)+ ?/ (12) 

Equilibrium is not achieved until L(t) reaches the size of the sample. In an infinite sample one thus observes an aging 
behaviour in the correlation function. Roughly speaking the system remembers its age t w through the value of its 
typical domain size L(t). 



A. The O(n) model 

Interestingly enough, one of the few exactly solved models of coarsening |l(J, namely the case of the 0(n) model 
with n large and a constraint 4> 2 = n, is also related to a problem which looks like a spin glass system. Indeed consider 
the following spin glass hamiltonian 

H = JijSiSj (13) 

ij 

where Sj are real spins with a spherical constraint sf — n, and Jij are random couplings. This model is usually 
called the (p = 2) spherical spin glass |^0| . Its Langevin dynamics, 

ds 

= E ■ 7 «*i(*) - Z W S ^) + »*(*) (14) 

3 

where z(t) is a Lagrange parameter inforcing the spherical constraint, can also be written in the basis where the Jy 
matrix is diagonal: 

^ = (X-z(t))s x (t)+r 1x (t) (15) 

Then the dynamical equation reduces to the one of the 0(n) model in Fourier space with A = — k 2 . The only 
important piece of information on the J matrix is the behaviour of its spectrum near its largest eigenvalue A* . The 
case of a square root singularity as for instance the Wigner law is equivalent to a d — 3 coarsening problem. Clearly 
the spherical spin glass does not really have a spin glass like behaviour (this has been known for long in the statics 
p0|); however it exhibits an aging dynamics, which is related to the growth of the correlation length in the O(n) 
model. Recently Cugliandolo and Dean |l5f| have performed a detailed computation of the dynamics of this problem, 
in which they computed all the relevant quantities of interest for our discussion. They found aging in the correlation 
function, but the anomaly of the response as defined in (1) vanishes, which shows that the memory effects are weak, 
besides, the overlap limit S(t w ) = lim^oo Q(t w +t,t w + 1) is a continuous function of t w , taking values between q 2 EA 
and qEA, with Soo = q E A, and, of course, lim tw ^ooQ(tm + t, t w + t 1 ) = C as (t + t 1 ). 

Physically, the system evolves in time through two processes: a thermal noise which affects evenly all the component 
s\, superimposed onto a deterministic evolution which amounts to reinforcing the eigenmodes closer to the A*, and 
this deterministic part is dominant over the thermal noise. 

The fact that this system has the property Soo — qEA together with the weak ergodicity breaking property 
limt^^oolimt^aoC^u, + t, t w ) = suggests the existence of a kind of gutter in phase space: two replicas, even 
decoupled, remember forever that they are evolving in the same canal. 



B. Scalar order parameter: analytic study 

We now turn to the domain growth problem in the case of a scalar order parameter. This problem cannot be 
solved exactly but we shall use a well known approximation [^8[ , recently developped by Bray and Humayun pT[ . 
We refer the reader to Bray's review jnj for a detailed presentation of the method. The idea is to take advantage of 
the universality of domain growth in the scaling regime: After an initial regime of fast growth, the order parameter 
saturates at the equilibrium value inside a domain and the only way for the system to further decrease the free energy 
is the reduction of the surface of walls between different domains. Therefore, the dynamical properties of the system 
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at late stage of growth are given by the motions of the walls and in particular by their curvature; the particular shape 
of the potential V[<f>], provided it has well separated minima, is not crucial. If the growth is influenced by an external 
field, the difference between the minima introduced by the field will be the relevant variable. The universality gives 
the freedom to choose an appropriate form for the potential in the free energy, and also a special form for the thermal 
noise, which make the analysis more tractable. Specifically the Langevin equation is replaced by: 



dcj>(x, t) 



where rj(x,t) is the gaussian white noise with zero mean and correlator: 

< r)(x, t)rj(x', t') >= 2T5(x - x')5(t - t') 
The field <f>{x,t) is parametrized by an auxiliary field m(x,t), through: 

1/2 



<j)[m] = (f>a 



dxexp — (x 2 /2) = 4>q erf[m/V2] 



erf' 1 








_4>a_ 





With the following choice of the two potentials: 



V Q [<t>] = ^exp-2 

7T 



the field m satisfies a very simple equation: 

dm 
~dt 



2 ,2 

, VM = ^expv^er/" 1 



V 2 m + (1 - (Vm) 2 )m + rj. 



(16) 



(17) 



(18) 



(19) 



(20) 



With the wall profile function ([Tq), the field m(x,t) measures the distance of the point x from the interface: at 
infinite distance from the wall, the field <f> saturates to its equilibrium value. Moreover, the potential Vq has the 
required two wells shape at the two equilibrium values. The choice of the potential V± does not alter this shape and, 
as can be seen from (|lq ) and from the fact that V{[<fi] = <fi' , it corresponds to a thermal noise acting only on the 
interface. This is an approximation which is not able, for instance, to reproduce a process of nucleation of a bubble. 
In other words, the value of qEA in this case remains fixed to the T = value, qEA = 4>o- However we expect that 
this approximation will not affect our main conclusions concerning the various large time limits of the overlap. 

The physical situation of a rapid quench will be represented by taking the boundary condition for m(x,t) to be 
gaussian with zero mean and correlator: 



< m(x, 0)m(x' , 0) >= 5(x — x'). 



(21) 



Equation ( p0| ) can be solved by neglecting the non linear term or, more correctly, by taking into account its mean 
value. Let us neglect it in a first approach. Equation (pp|) can then be solved, giving: 



m(x, t) 



\k\< 



d d k 
„-i 2tv*' 



ikx 



(22) 



where e is a cut off given by the width of the interface. The linearity of the equation and the independence of the 
boundary condition and the noise preserve the gaussian character of the probability distribution for the field to. Mean 
values of functions of the field cj> can be computed in term of the evolution of the first and second moment for the 
gaussian distribution of to. To compute the correlation function C(r+t w , t w ) we introduce two fields mi = m(x, T+t w ) 
and to 2 = m(x,t w ). When computing the overlap Q(t + t w ,T + t w ) the fields TOi and to 2 denote respectively 
to*- 1 -* (x, t + t w ) and to/ 2 ) (x, t + t w ). In both cases the joint distribution of mj and to 2 is a gaussian P(x±,X2), which 
we parametrize as: 



P(xi,x 2 ) 
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2wy/cr 1 (T2 

with <7, =< to 2 >, c 12 =< TO!TO 2 > and / = — | 
overlap) is given by: 



exp 



TL 
' 2 



x 2 
02 



2/a;ia;2 



(23) 



Vw 5 



V^2 

In terms of this distribution the correlation (or 



G 



r+oo 2 
< 0[toi]^[to 2 ] >= / dxidx2(p(xi)(f>(x2)P(x-i,X2) = 4>q — arcsin/. (24) 



where the function ( |Tq ) has been replaced by 4>[m] — 4>o sign[m], a good approximation in the large time regime. 
The calculation therefore reduces to the computation of the parameter / in the covariance matrix of the probability 
distribution P{x\,X2), which is easily obtained from (p2|). Defining : 



we obtain, for > 1: 



C(t +t w ,t w ) = O 



arcsin ■ 



1 + T^(2^,t + 2^) 



[1 + TT{2{t + ^), 2(r + t«,))]V2[l + T^(2i w , 2t w )]V2 j ( 26 ) 
and 

9 2 f 1 + TT{2t w ,2{r + t w )) ) , s 

The asymptotic behaviors of Q and C are very similar to those studied above in the domain growth of a n — ► oo 

component order parameter (or p = 2 spherical model) : the asympotic relation between Q and C at iu, — * oo with r 
finite is satisfied and, for fixed t w and r large, C has the limiting behavior 

■ x d/4 



r 



C(T + t«,t w )~ — (28) 



while Q does not go to zero and the limiting value S(t w ) is a continuous function of t w , approaching the equilibrium 
value qEA as t w grows: — cjea ■ So within this approximation this coarsening problem falls into the type I 
classification. 

Note: when one includes the effect of the gradient squared term in (^2|), treated as an average term (as in jn|), the 
result is similar except for a change in the numerical value of the function T in (p5|), where a term ccr~ 2 is present, 
instead of exp —a. 

In order to check this approximate analytic treatment, we have performed numerical simulations of domain growth 
in two dimensions, for a scalar field evolving with a Langevin equation, and also for Ising spins on a regular two- 
dimensional lattice, with Glauber dynamics pffl. 



C. Scalar order parameter: numerical studies 



We have simulated the evolution of a scalar field on a two-dimensional square lattice, according to the Langevin 
equation ([l2|), with a quartic Vq and a bold discretization scheme: 

4>{i, j,t+l) = <f>{i, j, t) + {cj)(i + l,j,t) + - 1, j, t) + <t>{i, j + 1, t) + (j>{i, j — l,t) 

-4 * cf>{i, j, t) + 4>(i, j, t) - <f>(i,j, t) 3 )*h + v (i, j, t) (29) 

where r\ is a gaussian noise with zero mean and variance 2Th, h being the used time-step. We proceed by parallel 
updating of the field, and vary the time-step h. At t = 0, 4>(i,j) are taken as independent random variables uniformly 
distributed between —1 and 1. We let the system evolve during t w according to (|2^), make a copy of it, and let the 
two copies evolve independently, i.e. with independent thermal noises. We record the correlation of each of the copies 
with the system at time t w and the overlap between the replicas. 

We present simulations at fixed temperature: we record the overlap and the correlation function for different values 
of t w , The linear size of the system was of 200 sites, and one run was made with a 400 * 400 lattice. Each simulation 
was made with three different values of h (h = .02, .04 and .08), to check that the results did not depend on the used 
time-step. We checked besides that the j- scaling is well obeyed for the correlation function for large enough t w ; for 
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the overlap, no such scaling is found. We plot the overlap at time t w + 1 versus the correlation between times t w and 
t w -\- 1. 

The Q versus C curves show quite clearly that the overlap, after a transient regime where it decays faster than the 
correlation, has a finite limit as C goes to zero. This limit grows with (3 and with t w . 




FIG. 2. Overlap Q(t w +t,t w +t) versus Correlation C(t w + t,t w ) for the scalar field in two dimensions, for (3 - 
waiting times (from bottom to top, 10, 20, 50 and 180 MC steps). 



6 and different 



This result agrees with the previous analytic study, as far as the asymptotic behaviour of the overlap and correlation 
are concerned. We have also performed simulations of a two-dimensionnal Ising spin system (with nearest-neighbour 
ferromagnetic interactions), with Metropolis dynamics with random updating R at each sweep through the lattice, 
spins are updated in random order, but this order is the same for both replicas. The results (see figure 3) agree with 
those obtained by Langevin dynamics. 



i * 



+ t * 

t t ; i * 



FIG. 3. Left:Overlap Q(t w + t,t w + i) versus Correlation C(t w + t,t m ) for the Ising model in two dimensions, for f3 — 2, 
t w = 2 3 (bottom) and t w = 2 4 (top); Right: Correlation and overlap for the Ising model in two dimensions, for different values 
of t w and t: t w = 2 3 , 2 4 , 2 9 , and t = 2, 2 10 

Note: We also made computer simulations for a model introduced in [ p5| , consisting of an Ising ferromagnet on 
a cubic lattice, with weak next-nearest-neighbour antiferromagnetic couplings; in this model, the growth is slowed 
from a power-law to a logarithmic behaviour; nevertheless, we find for the correlation and overlap functions a similar 
behaviour as for the simple ferromagnet. 



1 Notice that the choice of dynamics is important: as soon as the chosen algorithm is not deterministic at zero temperature 
(it is the case for example if we take Glauber dynamics), the overlap will decrease to zero even at T = 0. 
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D. The XY- model in one dimension 




A simple and soluble model where domain growth can be studied without approximations on the potential is the 
XY model in d = 1 pfj. Namely, the system has no phase transition for T > 0, but at very low temperature the 
correlation length L eq is very large. Then at time scales where the size of the domains is small compared to the 
correlation length, the system presents the typical non equilibrium features of a multiple phase system. 

In this simple model, the order parameter is a two dimensional vector field </)(x, t) of fixed length 4> 2 — 1 and the 
coarse grained free energy is: 

/ ^\ 2 

(30) 

Using the non linear mapping (j)(x,t) — (cos9(x,t),sui9(x,t)), the Langevin evolution equation for the field 9{x,t) 
can be easily written and explicitly solved without approximation. The physical situation of a rapid quench from 
a disordered phase to a very low temperature can be included in the formalism by taking the boundary condition 
9(x,0) to be gaussian with zero mean and correlated at distance £. As for the scalar order parameter model, the 
linearity of the Langevin equation preserves the gaussian character of the probability distribution for the field 9 and 
the problem can be solved by computing the time evolution of the moments. Let us now consider a quench to a very 
low temperature, i.e. a situation where the equilibrium correlation length L eq is very large, and the time needed for 
the domains to reach this size, t eq is also very large: in fact, L eq ~ T~?^ and t eq ~ T~ 2 . In a time regime where 
t + t w , t w 1 but L eq is very large compared to the size of the domains we have a very simple expression for the 
correlation function C(r + t Wl t w ) and for the overlap Q(t + t Wl r + t w ) of two replicas separated at time t w : 



C(t + t w ,t w ) = exp 1- [2(t + 2t w f/ 2 - [2(r + t w )Y' 2 - (2t w ) 1 / 2 

_ Il {2 [(r + 2t w y/ 2 - (rf' 2 ] - (2(r + t w )f' 2 - (2t w fl 2 } 



(31) 



and 



Q(r + t w ,T + t w ) =exp — 1 j=T{2t) 1/2 . (32) 

Since the size of the domain evolves as Lit) ~ ti, it is clear that there exists at very low temperatures a regime 
with 1 <C t w <C r -c t eq , where the correlation has already decayed to zero while the overlap has still a finite value. 
Indeed, C decays to zero with a term in the exponential that does not depend on the temperature, but the argument 
of the exponential for Q(t + t w , t + t w ) is L(t) 2 / L 2 q . 



E. Conclusion 



The previous study shows that the domain growth processes considered here are essentially deterministic in nature, 
and that their phase space is very simple: we indeed exhibit a time regime 1 <C t w <C t <gC t eq (t eq being the equilibration 
time, which is infinite in the true aging problems, but remains finite in the one-dimensional XY model) where the 
system at time t w + t has already drifted away from its position in phase space at t w (C(t w +t,t w ) is very small), while 
two copies separated at t w are still evolving together (Q(t w +t, t w + t) is finite). These type I systems are characterized 
by the existence of a finite limit Soo = lim^^oo lim^oo Q(t w + t, t w + t) (with S(t w ) = lim$_>oQ Q(t w +t,t w + t) 
growing continuously towards qea as t w grows). The system can therefore be thought of as moving along a gutter 
in phase space. It is reasonable to expect that, in these systems, the influence of the thermal noise will be limited 
in time, and there will be no anomaly in the response function (this has been shown so far only for the O(n) with 
n — * oo model pi). 



IV. TYPE II MODELS 
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A. A particle in a random potential 



We now turn to the study of a well known disordered mean- field model, where we expect a different kind of behaviour 
for the overlap: the toy model described by the hamiltonian p^,|5|: 

H = ^^2^ a + V(<h,.-,4>N), (33) 

a 

where V is a gaussian random potential with correlations: 



V($)V$) = -Nf P^fH (34) 



with: 



m = ~w^w (35) 

This model describes a particle in a random potential, in N dimensions, but it can also be interpreted as a spin 
glass model: the 4>a are then soft spins, in a quadratic well \^^ a and they interact via V; the statics has a low 
temperature spin glass phase, with continuous replica symmetry breaking for 7 < 1 (long-range correlations of the 
disorder) or one-step replica symmetry breaking for 7 > 1 (short-range correlations). Slightly different forms for / 
allow also to deal with the dynamical equations of the spherical p-spin |2^|I|,[32| model too. This system is described 
by the Hamiltonian ]33],[4| 

JV 



y ] Jii...i p s ii--- S i p (36) 



il<...<i p 



with the constraint s ? = 1> ano - gaussian distributed random p-spin interactions, ft can be described by a toy 

model, with 

/(&) — § C 1 -!)" ^ 



and a small modification of the dynamical equations (Al) and (A2), which amounts to implementing the spherical 
constraint by a time dependent Lagrange multiplier /z(£). 

-<1) -<2) 

To compute the overlap function, we introduce two replicas <p and 4> > with a Langevin dynamics: 



6 « (*) _ dH , .,«, 



st ~ ^V) '"^ (38) 

where 77 W and rj^ are two white noises with < r)£ (t)r) ( *) (t') >= 2T<W<Kt - i'), and r)£\t) = r/ ( i \t) if t < t w . For 
t > t w , rf 1 } and r/ 2 ) are uncorrelated. 

Using standard field theoretical techniques |l^,|3(]], it is now possible to derive the evolution equations for the 
correlation and response functions of each replica, C^'(t,t') = C^(t,t') = C(t,t'), and r^(t,t') = r^{t,t') = 
r{t,t'), and for the overlap Qt w (t,t'), in the large N limit. These quantities are defined by (||), (Q) and @ and the 
corresponding equations are written in appendix A. 

The large time limiting values of the correlation define q, qo, and q\ (sec p2[): 

lim^oo C(t, t) = q 

lim^oo C(t, t') = q Q (39) 
lim lim^oo C(t +T,t) = q\ 

r — >oo 

We will now study the behaviour of the overlap function in different time regimes. 

We first study the regime of asymptotic dynamics which corresponds to taking the limit t, if — ► 00, with r = t — if 
finite. We thus obtain the functions 
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s (r) = lim r{t' + T,t'), 

t' — >oo 



C as (r) = lim C(t'+r,t'), 



(40) 



? os (t,t') = lim Q tw {t w + T,t w + t). 



(41) 



In this regime, time-translational invariance and fluctuation dissipation theorem (FDT) are obeyed: Tr as {r) — 
~-jfpC as {T). It is well known that this asymptotic regime is identical to equilibrium dynamics for systems with long 
range correlations of the disorder 1 30 |||, but it is different for short range correlations ]29| , |30| ,p| . In both cases 
we have found, as expected from the general discussion of Sect. II, that Q as ( T i T ') = C os (r + rj. In particular, 
linv^oo Q as (T, t 1 ) = lirn r /_ (00 Q as {r, t') = q±. 

Let us now consider the aging regime. This regime corresponds to having time differences, like t — t w , diverge 
when t w — ► oo. Here we shall consider the overlap function Qt m {t,t') in the "double aging" regime where if — t w also 
diverges. 

There is no full solution of the aging regime in spin glass systems. What has been proposed so far, in all cases, is 
an Ansatz about the behaviour of the correlation or response. The first such proposal, by Cugliandolo and Kurchan 
H , concerns the case of the p-spin model. They showed that the dynamical equations can be solved in the long time 
regime (where one can neglect the time derivatives in ( Al) ) by the Ansatz: C(t, if) — C{if /t), r(t, t') = {x/T)C {if /t), 
or actually by any solution obtained from this one through a reparametrization of time t — > h{t), with h an arbitrary 
increasing function. This solution was subsequently extended to more complicated problems in which the static 
solution involves a full rsb, like the toy model with long range correlations of the noise jg) and the SK model Q. The 
case of the toy model with short range correlations of the noise has also been studied recently |32| . The formalisms 
developped in j^] (non overlapping time domains) and in || (triangular relations) represent the same Ansatz but look 
rather different. Here we shall present the Ansatz using mainly the former approach, together with the necessary 
ingredients for understanding the correspondence between the two formalisms. 

Considering first two time quantities like the correlation or response, the aging regime corresponds to sending i and 
tf both to infinity, the difference t — if being itself divergent in this limit. The dynamical equations can be solved in 
this limit (up to a time reparametrization), neglecting the time derivatives. We consider non overlapping domains 
in the {t, if) plane: two times t and i', with if < t, belong to the same domain D u if we take the limits t — > oo, 
if — ► oo, with the ratio h u {t')/h u {t) finite and fixed to e~ T (0 < r < oo). The h u are a family of increasing functions 
indexed by a parameter u, such that, if w < u < v and the times t, t' belong to D u , then h v {t')/h v {t) = and 

h w {t')/h w {t) = 1 (a possible choice is h u {t) = exp{t u ), in which case D u is such that if = t — r). The domain 

U = 1, with hi{t) = exp(t), corresponds to the asymptotic regime where FDT and TTI hold. We find convenient to 
express the fact that (t, if) G D u by the following diagram: 




FIG. 4. {t,t') £ D u , with h u (t')/h u {t) = e 



It is then easy to show that, if we consider three times if < s <t, with (s, t) G D u and {if, s) G D v , then {if, t) belong 
to D rain („ „) (see figure 5) which is an ultrametric inequality. If for example v > u, we have indeed h u {t')/h u {s) = 1, 
so h u {t')/h u {t) = h u {s)/h u {t). 
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In each domain D u , we assume the correlation and response to behave as: 

C(M') = C u (r), r(t,t') = dHh ^ ]) T u {T), (42) 

with a continuity condition: if D u a nd D v are neighbouring domains, with u < v, then C u (0) = C v (oo). Then it 
is possible to rewrite the equations (Al) (see appendix B for details), and to show |3lJ] that they possess solutions 
obeying a generalized form of the fdt relation called "quasi-fdt" : 

Xu ^pL = _Tr„(r) (43) 

If we now consider three times, t' < s < t (see figure 5) the Ansatz implies a simple relation between correlations 
at times t, t', t, s and if, s. When (if , s) and (s, t) belong to two different domains we have: 

C(t, t') = min(C(t, s), C(s, t')), (44) 

and if they are in the same domain D u , with 

^ = e -,M£i = ^, tl „MG^— , (45, 
h u (s) h u (t) h u (t) 

so that 



C(t,t') = C u (r' + T)=j- 1 (j u (C(t,s))j u (C(s,t'))), (46) 

where C(t, s) — C u (t), C(s, if) = C u (t'), and j u (z) = cxp(C^ 1 (z)). Equations ( ^4] , p6| ) form the basis of the formalism 
of triangular relations introduced in J(| , and applied to the toy model in |}2| . We provide in appendix C the solution 
for the overlap function using this formalism. 

Since the overlap function Qt w (t,tf) involves three times, we are now looking for a function depending on the 
domains D u and D u >, Q„, tt '(r,r'), where (t w ,t) <E D u and (t w ,t') € D u > {h u (t w )/h u (t) = e _r , h u > (t w ) / 'h u > (*') = e~ T '). 

We rewrite in appendix B the equations (A2) in this frame, and show that they are solved by the following Ansatz: 
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If u ^ u' (for example u < u'): Qt w (t,if) — C u (t) — min(C(t,t w ),C(t' ,t w )) . 

K« = W: Qtjt,t') = C u (r' +t) = jZ X {ju(C{t,t w ))j u {C{t' \t w ))) . (47) 

This Ansatz can be easily understood in terms of the previously introduced diagrams: at t w two "time-sheets" 
separate (see figure 6) and two ultrametric systems appear, one for each replica. We stress that this solution exists 
independently of the actual choice of the disorder correlation, and therefore it is independent on the precise solution 
of the aging dynamics: whatever the number of non overlapping domains appearing in this solution, whatever the 
actual solutions C u (r), there exists a solution for the overlap function in the aging regime which is related to the 
correlation by ([47|). 




Depending on the model, the variable u can be a priori continuous, or discrete. It was shown in |J, using the results 
from fjlfl , that for the long-range model u becomes a continuous variable. For short-range models on the contrary, 
which exhibit statically a one-step replica symmetry breaking, it has been shown [p[|3^] that the aging dynamics is 
solved by using a single time domain D u * (beside the fdt domain D\). 

To summarize, we have shown that, in the aging regime: 

• for the long range model: 

Q tw (t,t') = mm(C(t,t w ),C(t',t w )) (48) 

In particular, it is then clear that the long-time limit of Qt w {t,t') is qo] 

• for the short-range model, there exists a function j such that 

QtJt,t')=r 1 (j(C(t,t w ))j(C(t',t w ))). (49) 



Since j(qo) — 0, and since lim^oo C(t, t w ) = limt'_, 00 C(t' , t w ) = qo, we also have limine Qt w (£, t') = 
linv^oo Q tw {t,if) = q 

In both cases, limt^—^ limt_, 00 Q tw (t, t) is qo, different from lim^oo limt^^oo Qt w (t,t). Besides, no finite waiting 
time is sufficient to give a higher limit than qo for lim^oo Qt„(i,t): an increase in the waiting time only slows down 
the dynamics, but has no effect on the limiting values. No continuous approach to equilibrium can thus be seen in 
this way. 

Note: since the equations ( Al A2) are causal, a numerical integration is available, as in ||. Nevertheless, the integro- 



diffcrential character of these equations makes it difficult to reach very long times (a huge amount of computer memory 
is needed). Therefore, the numerical integration we were able to realize, although fully compatible with the previous 
study, were not conclusive enough to confirm it. 

For the case of the p-spin spherical spin-glass model, an analytic solution of the equations is also available: they 
are solved by the Ansatz corresponding to the short-range model, equation (|^). Indeed, we propose for the aging 
regime the Ansatz Q tul (t, t') — Q(t w /t, t w /t'), and we find that the equation giving Q tw (t, t') can be rewritten so that 
the three times t w , t and t' appear only through the ratios t w /t and t w /t', and that Q tvj {t,t') — ^C(t,t w )C(t' ,t w ) is 
solution of this equation. 
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We have thus shown that the overlap in the p-spin (p > 3) spherical model exhibits aging in a similar fashion as 
the correlation function, and decays to zero for any finite t w (qq = for this model). This behaviour is thus very 
different from the p = 2 or the domain-growth case. 

B. Aging in traps 

A trap model was introduced in Jl6| and developped in Jl7| to reproduce off equilibrium dynamics in glassy systems, 
and aging. The model consists of ./V traps with exponentially distributed energy barriers. This distribution leads to 
trapping time with infinite mean, and thus to aging. 

In the simplest version, Jl6[ ] the basic object is IIi(i, t w ), probability that the system has not jumped out of his trap 
between t w and t w + t. The overlap between two different states is zero, and the self-overlap is Qea- The correlation 
function is then qEA^-i(t,t w ). 

We now deal with two systems after t w : we introduce 11^ (t, t ,t w ), probability that the first replica has not jumped 
between t w and t w + t, and that the second one has not jumped between t w and t w +t' . The overlap Qt w (t w +i, t w + t') 
is then simply qEA^i [t,t ,t w ). 

If the system is in trap (3 (of lifetime 173) at time t w , this probability is e T " e t p , and thus we get 11^ (t, t , t w ) = 
II 1 (t + t',t w ). The overlap between the replicas is therefore: 

Q tui (t w + t,t w + t r ) = C(t w +t + t',t w ) (50) 

If we introduce a multi-layer tree, the only difference is that we now have a set of IL, (t , t w ) (j = 1 , . . . , M) , probability 
that the system has not jumped beyond the j th level of the tree between t w and t w + t. It is then clear that the 
equation (|50| ) is not changed, although the analytic expression for the correlation function depends on the parameters 
of the tree. 

For this particular model, the equilibrium relation is in fact satisfied even for out of equilibrium dynamics, because 
of the properties of the chosen exponential decay from the traps. It is then clear that the overlap function goes to 
zero for large times, for any finite t w , since lim^oo C(t w + 1, t w ) — 0: S(t w ) — 0. We have thus the same scenario as 
for the toy-model: for any finite t w , the overlap decays to its minimum allowed value, while an infinite t w gives qea 
as a limit. 



V. CONCLUSIONS 

In this paper, we have shown that the overlap between two copies of a system, identical until a waiting time 
t w , and then totally independent, is a quantity of interest regarding the geometry of phase space. We have indeed 
studied this quantity for several models, and shown that its decay is intimately related to the complexity of the 
landscape and to the type of aging. For simple systems, the long time limit of the overlap can be put closer and 
closer to the equilibrium limit qEA by changing the time the replicas spend together. On the contrary, for systems 
exhibiting a complex phase space, the limit of the overlap is always the minimum value, i.e. the two replicas are 
able to separate from each other, no matter how long the waiting time is (if it stays finite). This difference of 
behaviour can be quantitatively seen by a study of the various long time limits of the overlap: for any system, 
lim^oo lim^-Kx, Q tu , (t w +t,t w + 1) is qEA, but the inverse order of limits, Soo = liirit^oo lim^oo Q tw (t w + 1, t w + 1) 
(and the behaviour of S{t w ) — lim^oo Qt w {t w +t,t w + t)), distinguishes between aging in a "simple" phase space, 
and type II (spin glass) systems. 
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APPENDIX A: DYNAMICAL EQUATIONS FOR THE TOY-MODEL 

In the limit of infinite N, the dynamical equations for the two-times correlation and response functions (with t > t') 
read g: 
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dr(t t') f* 

g t ' = -l*r(t, f) + / ds m(t,s)(r(t,l/) -r(s, f)) , 

dC(t,t') 



at 



1 dC{t,t) 

2 ~dt 



and the equations we obtain for Q tw are: 



-fj.C(t, t') + 2 [ ds w(t, s) r(t', s) 
Jo 

t 

ds m(t,s) {C{t,t') - C{s,t')) , 

o 

-flC(t,t)+2 [ dsw(t,s) r(t,s) 
Jo 

t 

ds m{t 7 s) {C(t,t)-C{s,t)) +T , (Al) 



dQtJt,t<) = _^ (M , ) + 2 f l dsWt j tiS) r{t ' iS) 



i) 



at 

+ f dsm(t,s) (Q tw (t,t') - Qt.(*,0) 



o 

t 



with the same notation as [El : 



and: 



o 

+ / dsm(M) {Qt w (t,t)-Q tw (s,t))+Te(t w -t) , (A2) 
Jo 



«;(*, t') = f(b(t, 0), m(t, = 4/"(&(t, f )M«, 

6(t, t') = C(i, t) + C(i', t') - 2C(t, t') (A3) 



W tw (t,t') = f(B tw (t,t')) 

B tw (t, t') = C(t, t) + C(t', t') - 2Q tw (t, t') . (A4) 



APPENDIX B: NON OVERLAPPING TIME DOMAINS 



We show how to compute one of the integrals of equations (Al,A2) in the frame of the Ansatz described in [gj, 



using the diagrams introduced in figure 4. We apply this method to the overlap equation and show that a similar 
Ansatz is solution. 

For (i',t) S D u , we parametrize 

b{t,t>) = b u (r), r(t,t') = dHk d f )] r u {r), 
. d\n(h u (t')) 

m(t, t ) = — m u (T), w(t, t ) = w u {t). (til) 



The integral J., ds m(t, s) r(s,tf), appearing in (Al), has then three contributions (see figure 7): 



• (t', s) S D v (v > u), where m(t, s) = m u (r) d ln ^ s " — , 

• t' t, s in the same domain D,. 



'U ! 



(s, t) € D v (v >u), where r(s,t') = r u (r). 
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FIG. 7. For t' < s < t, three regimes can be separated. 

f* . , n / a dlnhJt) f . ^ f°° , , , r , , , , 
/ ds m(t,s) r(s,t ) = — | m u (T) } j I da r v (aj + / da m u {a) r„(r - a) 

Jf M V v>u J JO 

+ r «( T ) da m v (a) 

v>u J J 

Separating in this way all different contributions in the integrals, we obtain the following equation for b u (r): 



(B2) 



= K(t) 



/>00 PT 

M + JZ / ds m v (s) + 2T- / ds m u (s) b u (r - s) 

v<u J " j J " 

-4w u (t) y~] / ds r v (s) - / ds [m u (r + s) b u (s) + 4 w u (t + s) r„(s)] 

P OO 

+ ^2 ds I TO «( S ) M s ) + 4 w v(s) r v (s)} , 



(B3) 



The same method can be applied to the equation for the overlap function B tul (t,t') = B UiU >(t,t'). For u, v! < 1 
(with u ^ u') it reads: 



= S„,„/(t ) r') 



-(J. + ^ / ds m„(s) 

OO 



+ 2T- [ dsm u (s) B UiU ,(t-s,t') 
Jo 

r°° _ /*°° 

-4io u (t) ^ / ds r„(s) -4 W tt , u '(r,r / ) ^ / ds r v (s) 

u'>v>v,J° v>u'J° 

/OO /■T 
ds rv (s) - 4 y ds W U)U ' (r, r' - s) r„/ (s) 

/•OO 

ds [m u (r + s) 6 u (s) + 4 to„(r + s) r u (s) - m„(s) 6„(s) - 4 io u (s) r„(s)] 

Jo 
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/•OO 

^ / ds [m v (s) b v (s) + 4 w v (s) r v (s) 



If we insert in this equation the Ansatz 

B u , u >(t,t') = b u {r) (i{u<v!) 



we reobtain equation (B3). For u = u , the equation is slightly different, and the Ansatz 

B U , U (T,T') =b u {T + T'), 

together with the quasi- fdt, again gives back equation (B3), evaluated at r + t' . 



(B4) 

(B5) 
(B6) 



APPENDIX C: TRIANGULAR RELATIONS 



We derive the solution for the overlap function using the for ma lism of triangular relations ||,|3^]. Neglecting time 
derivatives, and distinguishing aging and FDT regimes in (A1,A2), the final equation for the slow varying part (i.e. 
non fdt) of B tw (t,t') = C(t,t) + C(f ,tf) - 2Q tw (t,t?) reads (with q = 2(q - qi j): 

2q 



fJ>+ ds m(t, s) 
Jo 

t rtw ft' 



IT - Y lW t Jt,t')-f(q)} 



= B tw (t,t') 

+ 4/ ds w(t, s) r(t, s) - 4 / ds w(t, s) r(t' , s) - 4 / ds W tvi (t,s) r(t',s) 



+ / ds m(t, s) b(t, s) — / ds m(t, s) b(t? , s) — / ds m(t, s) Bt w (t' , s) 



(CI) 



where now all the times in the equations belong to the aging regime. 

The approach of triangular relations measures the various time domains directly in terms of the distance b(t,t'). 
Indeed, there is a one to one correspondence between b(t,t') (at large times) and the functions 6 u (r). More precisely, 
if t,t' — ► 00 with b(t,t') fixed to B, then t,t' belong to D u , with h u (t')/h u (t) = e~ T , where u and r are fixed by 
bu{r) = B. Then, the quasi-fdt can be expressed as 



r(t,t')=X[b(t,t')}^p- 



(C2) 



where X[b(t,t')] — x u /(2T). Similarly, the ultrametric structure of time domains described in ( ^4] , ^6| ) can be written 
in a compact form as a triangular relation [^],[32| : 

b(t,t')=g(b(t,s),b(s,t')), (C3) 

with g(b,b') — max(6, b') when b and b' belong to different domains (also named "blobs" ||), and g(b,b') — 
3~ (i(fr)i(fr')) within the same domain. 

The Ansatz concerning the function B tm reads: 

B tw (t,t') = j[b(t,t w ),b(t',t w )} (C4) 

where 7 is a function to be determined. Setting 6(t, t w ) = b w , b(t' , t w ) = b' w , 6(t, t') — b, bo = 2(q — qo) and 



F(b) = - f ds X(s) 

Jb 



(C5) 



we obtain: 



= j(b w ,b'J 



+ 2T-^[/'[7(^,y-/'(?)] 



-M + 4 / dsf"(s)X(s) 

J b 

+ 4/'(&o) F(b ) +4 f ds f"{s) X(s) s + 4 f ds f(s) X(s) 

J b J bo 
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ds f"(s) F[g(b, s)]-i ds f"(s) X(s) g(b, s) 

J b 

4 ^ ds F\g(b,s)] f"[^(b w ,g(s,b w ))W(b w ,g(s,b w ))g'(s,b w ) 

4 j" ds f"{s)X{s) 1 [g{s,b w ),b' w ] (C6) 



where g is the reciprocal function of g: given three times t' < s < t, in the limit t' , s, t — > oo, 

b(t, = ff (6(t, fl ), 6( a , 0), &(*, = §(&(*, s),b(t, t')). (C7) 

Equation ( |C6| ) is a functional equation that gives 7 in terms of /. X and 5. 

Using the g-function described above, one can check from (|Cq ) that the solution is 7(x, y) = g{x,y) (this means 
that the relation between the overlap function Bt w (t, i') and the correlation functions r^,) and b(t' , t w ) is the same 
as the triangular relation between the correlation functions) . This way, one gets back the result for the overlap given 
in pp|). 
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